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SPHERICAL BALLOON RESPONSE TO THREE-DIMENSIONAL 
TIME-DEPENDENT FLOWS 

SUMMARY 


To fully understand what a balloon wind observation means, it is 
necessary to know the response properties of the balloon relative to the im- 
posed wind field. In this report, the basic equations of motion for a spheri- 
cal balloon are derived with the assumption that the atmosphere affects the 
balloon, but the balloon does not affect the atmosphere. The response of the 
balloon to the wind field is affected by aerodynamic drag and lift forces, vir- 
tual mass effects, the Archimedean buoyancy forces, and dynamic buoyancy 
forces. An analysis shows that the Coriolis terms in the equations of motion 
can be neglected for most rising or falling balloons. 

The equations of motion for a spherical balloon are linearized about 
a terminal velocity state in which the horizontal velocity is in equilibrium 
with the mean horizontal velocity of the atmosphere, and the rise or fall rate 
is that which is obtained when the Archimedean buoyancy force balances the 
aerodynamic drag force. The mean atmospheric flow is assumed to be a 
space -time invariant horizontal flow. The atmospheric flow and the balloon 
velocity vector are perturbed by superimposing small perturbations such that 
the square and higher-order terms in perturbation quantities can be neglected 
with respect to first-order terms. The environmental perturbation velocity 
field is a three-dimensional time -dependent vector field, which is represented 
as a four-fold Fourier integral involving three mutually orthogonal wave 
numbers and a frequency. The components of balloon velocity are, in turn, 
represented in terms of Fourier integrals which involve a frequency that is a 
function of the wind field wave numbers, frequency, and the velocity vector 
of the unperturbed wind field. The combination of these Fourier representa- 
tions with the equations of motion results in expressions which relate the 
Fourier amplitudes of the balloon to the Fourier amplitudes of the wind field 
via a three-fold integral. This integral convolves the wind field Fourier am- 
plitudes with the ballon system function over the domains of the wind field 
horizontal wave numbers and frequency. 



The nonstationary spectra of the balloon velocity vector are calculated 
in terms of the nonhomogeneous and nonstationary spectra of the wind field. 
The results are simplified by assuming that the wind field is statistically 
stationary and homogeneous. The stationary balloon spectra are thus obtained 
as three-fold integrals of the wind field spectra modified by the balloon trans- 
fer function. The introduction of Taylor’s hypothesis reduces these integrals 
to algebraic expressions which show that the spectra of the corresponding 
components of velocity of the balloon and wind field are proportional and that 
the function of proportionality is a function of the balloon parameters and the 
frequency of the wind perturbations relative to the balloon. 

INTRODUCTION 


It is widely recognized that meteorological balloons do not respond 
perfectly to the wind fields they traverse. A number of analyses have been 
published during the last 10 years which attempt to define the nature of balloon 
responses to the wind. Among these works are those by Reed [ 1] , Lewis and 
Engler [2] , Zartarian and Thompson [3] , Eckstrom [4] , Scoggins [5] , and 
Fichtl [6] . However, in all cases the analyses have an incomplete set of 
equations of motion as their starting point. In addition, the analyses are 
restricted to rather simplified atmospheric flows in which only vertical varia- 
tions of the wind are permitted. This report presents what are believed to be 
complete equations of motion for a spherical balloon and also presents a 
linear theory of spherical balloon response to three-dimensional, time- 
dependent atmospheric flows. 

EQUATION OF MOTION 


The equations of motion for a spherical balloon of mass m (skin and 
inflation gas) relative to a Cartesian frame of reference attached to the earth 
with the x 3 -axis directed toward the zenith and the x A and x 2 axes in a plane 
tangent to the earth (e.g. , directed toward the east and north) are given in 
tensor form by 


/ dv k 

m V "dT + 2!2 i V j 6 ijk/ 


“• gm{ k3 tD k tL k 
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+ m 


3T ( \'V +c i (u i' v i ,e ii 


ijk 
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^ ’ 


( 1 ) 


where t denotes time and all subscripts and superscripts take on values 1, 

2, 3. We shall use the Einstein summation convention; i.e. , repeated sub- 
or superscripts imply summation from i to 3. The quantities v^ and 

note the k components of velocity of the center of gravity (eg) of the balloon 
and the atmosphere at the location of the eg in the absence of the balloon. 



The first term on the right-hand side of this equation denotes the 
gravitational body force on the balloon, where g is the magnitude of the 
gravity vector and 

1, k = 3 

0, k ^ 3 . 

The second and third terms denote the k components of the aerodynamic 
drag and lift forces. The drag and lift forces act parallel and normal to the 
relative wind vector u^ - v^. The standard assumed forms of these forces 

which have been used in analyses of balloons are 

D k = k p AC D ^ Ul " Vl ^ 2 + * U2 “ V2 ^ 2 + * Us " V3 ^ ^ (\“\) » 

L k = 2 p A C Lk [ ( Ul " Vl ) 2+ K- v 2) 2+ ( u 3 “ v 3) ? l (lb) 

where p is the density of air, is the aerodynamic drag coefficient, and 
C Lk iS the k component of the aerodynamic lift coefficient vector. All avail- 
able evidence appears to show that these representations are valid for the 
Reynolds number region of typical balloons; i.e., iO 4 < Re < 10 6 , where 
Re = v 3 D/y , v 3 being the mean rise or fall rate, D the balloon diameter, 
and v the coefficient of kinematic viscosity. 

The fourth term is the force which results from the apparent or vir- 
tual mass effect, where m a is the apparent mass and O is the i component 
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of the angular velocity vector associated with the rotation of the earth. The 

quantity e is the alternating tensor which has the following properties: 
ijk 


e 


ijk 


1 , i ^ j ^ k, even permutation 
-1 , i f j f- k, odd permutation 

0 , otherwise. 


The fifth term is the force on the balloon resulting from atmospheric 
pressure gradients, 9p/8x^, where V is the volume of the balloon. It is 

reasonable to hypothesize that the pressure gradient force on the balloon 
should be proportional to the atmospheric pressure gradient at the balloon eg 
and the volume of the balloon. However, it should be remembered that the 
pressure gradient is that of the atmophere in the absence of the balloon. As 
in the specification of the atmospheric velocity vector at the balloon eg, it is 
assumed that the atmosphere affects the balloon but the balloon does not af- 
fect the atmosphere. Thus, the form of the pressure gradient force in equa- 
tion (1) is a hypothesis. However, as we shall see, it leads to the accepted 
form of the Archimedean buoyancy force for balloons. 


The net Vector sum of the forces on the right side of equation ( 1 ) must 
be balanced by the inertial reaction, which is the left side of this equation. 
The two Coriolis terms in equation (l) which contain £2 result from the fact 

that the frame of reference, attached to a point on the surface of the earth, is 
not an inertial frame. These terms are a natural consequence of the trans- 
formation of a Lagrangian time derivative of a vector relative to an absolute 
or inertial frame of reference to an accelerating frame of reference. 


ELIMINATION OF THE PRESSURE GRADIENT FORCE 


The inviscid equations of motion of the atmosphere will be used to ob- 
tain an expression for the pressure gradient 8 p/9 x^. These equations are 
given by 


Dt 


1 

P 


d£ 

9x, 


- 2 


£2 . u. e... 

i j ijk 
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3 k 


( 2 ) 
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Elimination of the pressure gradient between equations (i) and (2) yields 


(m + m ) 
a 


dv 

— + (2m + m a ) S2.VJ e [jk = - g <m - mj 6^ 


+ D. + L, + m 
k k a 



S2 . u. 
i J 



+ m 
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D \ 


Dt + 2!J i U j e 
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(3) 


where m Q (= pV) is the mass of the atmosphere displaced by the balloon. The 

first term on the right side of this equation is the well-known Archimedean 
buoyancy force which results from combination of the gravitational body forces 
in equation (i) and (2) . The last term, the dynamic buoyancy force, results 
from inertial accelerations of the atmosphere. 


At this point, a brief discussion of the various derivatives in equation 
(3) is appropriate. The derivative in the last term on the right side of equa- 
tion ( 3) is the material derivative following the atmospheric particle at the 
balloon eg. It can be expressed as 


Dt 


^ 


at 


+ u. 


ax. 

j 


(4) 


To calculate Du^/Dt, we perform the Eulerian operations on the right side of 

equation ( 4) and then evaluate at the location of the balloon eg. The differenti- 
ation indicated on the left side of equation (3) and in the apparent mass term 
on the right side denotes a change following the balloon. To evaluate the de- 
rivative of u. in the apparent mass term, we note that 


^\ = ^k 

dt 9t + V j 


*k 

9x. 

J 


(5) 


Thus, to determine the change of u^ following the balloon, we execute the 
Eulerian operations on the right side of equation (5) and then evaluate at the 
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location of the balloon. The difference between the operators d/dt and D/Dt 
becomes obvious by comparing the right sides of equations ( 4) and ( 5) . 

Utilization of equations (4) and (5) permits equation (3) to be expres- 
sed as 


dv k 

(m + m a> "dT + (2m + m a } fi i V j e ijk = -S ( m - m 0 } 6 k3 


d \ ^k 

+ D, + L, + (m + m ) — — — + (m v + m u ) — — 
k k a o 0t a j o i 9x. 

J 


+ (m + 2 m ) n. u. e. 
a o l j ljk 


(6) 


GENERALIZED RESPONSE PROBLEM 


The motivation behind the formulation of the equations of motion for a 
balloon is the calculation of the responses of balloons to atmospheric fields 
of velocity. Accordingly, this permits one to determine how well a balloon 
can measure the wind and what that measurement means. We denote the 
Lagrangian displacements of the balloon eg relative to the frame of reference 
attached to the earth with X (t) . The Lagrangian displacements are functions 

K 

of time and by definition are related to the balloon components of velocity 
through the expression 


V*» - 



(7) 


Substitution of equation ( 7) into ( 6) yields 
d 2 X^ dX. 

(m+m a ) — + < 2m+m a > a tl* e .i)k = -S (m - m „ )s k3 




dX 


m 


9u, 


j. , u I + (m + m ) — - 
dt m / a o dt 
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dX. 

-jJ- + (m a + 2m o ) n, Uj e ljk . 

( 8 ) 

This tensor equation represents three equations for the balloon coordinates. 

To solve these equations, we must specify the initial conditions and the 
wind field. Substitution of the known wind field into equation ( 8) and evalua- 
tion at x = X (t) yields the equations for the Lagrangian displacements, 
m m 

The solution of these equations yields the functions X m (t) and hence 

(v * dX /dt). The balloon velocities v (t) can be compared with the speci- 
v m m m 

fied winds along the trajectory u^ [X^ (t), t] to determine how the balloon 
responds to the prescribed wind field. 


+ m u. 
o J 


ax. 

] 


+ m 


ax. 

3 


The inverse problem is to specify the functions X (t) for a single bal- 

m 

loon track and determine the wind field (t) , t] along the balloon track. 

Equation (8) is a partial differential equation in u^, and because it is a partial 

differential equation, the stated inverse problem cannot be solved. To solve 
this equation for u^, one must do the impossible; namely, release an infinite 

number of balloons so that can be determined as a function of time and 

space. One could then solve equation ( 8) for u m ( x j ? » t) and then evaluate for 

the particular balloon in question to obtain u [X^(t) ,t] . 


SOME SIMPLE BALLOON PROBLEMS - DEFINITION OF A 
BALLOON TIME CONSTANT 


Let us consider the simple situation of a vertically rising or felling 
balloon in a quiescent nonrotating atmosphere. Furthermore, the aerody- 
namic lift is assumed to be zero. Under these assumptions, equation (6) 
reduces to 


dV 3 i 

(m + m a ) — jj- = - j P AC d |v 3 I v 3 + g (m Q - m) 


Let us assume 


(9) 


v 3 = v 3 + vj 


( 10 ). 
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where v 3 is a mean rise or fall rate in equilibrium with the mean flow drag 
force and buoyancy force, such that 

PAC D |v 3 | v 3 = g (m Q -m) . (11) 

The quantity v 3 is a superimposed perturbation which is sufficiently small in 
magnitude such that second-order terms in v 3 can be neglected with respect to 
first-order terms. Substitution of equation (10) and (11) into equation (9) 
yields 


(m + m a > = -|p A c p Jv, 1 v, . (12) 

Thus, we are considering small balloon velocity perturbations superimposed 
on a mean terminal vertical rise or fall rate. Integration of equation ( 12) 
yields 


Vg (t) = V ; (0)e" t/T , 

where 

2(m + ro & ) 

T= pAC D lv 3 l ’ 

or in view of equation (11) 

m + m |v | 

a 3 

|m-m | g 


(13) 


(14a) 


(14b) 


The quantity T is the balloon time constant and is that time required for an 
initial balloon vertical velocity perturbation from terminal velocity conditions 
to be reduced by a factor of e -1 . The balloon time constant is in the order of 
a few seconds for most meteorological balloons, and it commonly occurs in 
various theories of the responses of balloons to atmospheric motions [6, 7] . 

For a neutrally buoyant balloon, the time constant is defined in terms 
of the solution of equation ( 9) for m = m Q . We take v 3 to be initially positive 
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[v 3 (0) > 0] and assume that v 3 is positive definite, so that upon being dis- 
placed the balloon will immediately decelerate as it passes through the atmos- 
phere. This means we can drop the absolute value sign in equation (9) and 
obtain the solution 

v 3 (t) i 

v 3 (0) i+t/T ' 1 1 

When t/T = 1, we find from equation (15) that v 3 (t)/v 3 (0) = 1/2. Thus, the 
quantity T is the time required for the initial velocity of a neutrally buoyant 
balloon to be reduced by a factor of one -half when the only effect present is 
that of aerodynamic drag. 


THE CORIOLIS TERMS 


Let us consider the Coriolis terms in equation (6) . ’to neglect the 
Coriolis terms on the left side of equation ( 6) we could require 


dv i 

dt 

V 3 - fi 3 v 2 | 


> > 


2m + m 


m + m 


(16) 


a 

to 

2m + m 

I dt 

a 


IfisVt! >> m + m a 


(17) 


g m - m 2m + m 

0 a , . 

1 1 TT^ r >> ; • (18) 

(m + m)|fi 2 v 1 | m + m 

a. cL 

In the analysis that follows, we shall set (2m + m )/(m + m ) equal to its 

SL cL 

maximum value of 2. If these inequalities are satisfied for this value of 
(2m + m )/(m + m ) , they will be satisfied for all values of this parameter. 

cl cl 

In the case of an ascending or descending balloon, the vertical velocity should 
be near terminal conditions, so that we can cast equation (18) into the form 


9 


T fi 2 < < 


(19) 


2 'V 

Now, |v 3 1/ !vj | is in the order of unity, T has value approximately equal to a 
few seconds, and J2 2 ^ 10“ 4 s -1 . Thus, condition (19) is satisfied for most 
rising or falling balloons. 

To estimate the left side of equation ( 17) , we use the frequency o 

B 

of the wind perturbation relative to a rising or falling balloon to estimate the 
operator d/dt in equation (17) and thus require 

fi B 2lv i* 

1071'"“ • (20 > 

Now, |v 1 |/v 2 l~ 1 so that equation ( 20) requires 12 B /|n 3 [ » 1. The quantity 

ft is on the order of 2 tt x 10“ 2 for a balloon with rise or fall rate equal to 
B 

approximately 5 m- s -1 ascending or descending through a wind perturbation 
with a wavelength approximately equal to 500 m and |ft 3 | ~ 10~ 4 s“* , so that 
ft g/ ft 3 2ir x 10 +2 . Thus condition (20) or (17) appears to be satisfied for 

typical perturbations on the wind profile. 

Let uS%ow consider condition ( 16) . For a rising or falling balloon, 
v 3 tends to have one sign with magnitude approximately equal to the terminal 
velocity. It is highly unlikely that ft 2 v 3 = ft 3 v 2 for all t, so that the possi- 
bility of the denominator vanishing in condition (16) will only occur at a few 
brief instants at most. If the denominator vanished, the x 1 component of the 
Coriolis effect would then be zero, and we would not have to be concerned 
with condition ( 16) . Since 


I2 2 |V 3 I + I n 3 v 2 l S: ln 2 v 3 -fi 3 v 2 | 
we replace condition ( 16) with 


fi 2 lv 3 l + l£2 3 v 2 l 


( 21 ) 
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where we have used n^lVjl as an estimate of IdVj/dt I . Because fi 2 — | £ 2 3 l 
in mid-latitudes and IvJ — lv 2 l, we have 


0 , 


* 0,1 


>> 2 



( 22 ) 


It is clear from the previous comments that this expression is satisfied for 
the rising or falling balloon. 


Let us now consider the Coriolis terms on the right side of equation 
( 6 ) . To neglect these terms we could require that 


dv. 


dt 


I SI 2 u 3 ~ «3 U 2 | 
dv„ 


dt 


Uil 


> > 2 


> > 2 , 


(23) 


(24) 


and 


g - m | 
o 

2 (m + m) £22 I | 

cL 


> > 2 


(25) 


Since luj ~ Ivjl, conditions (24) and (25) are satisfied because the 
replacement of lujf with (v t | yields conditions (17 ) and ( 18) , which have 
already been shown to be satisfied by the balloon motions. We replace con- 
dition ( 23) with 


£2 


B 1 

fi 2 I u 3 1 + 1 £2 3 u 2 1 


> > 2 


(26) 

9 


as in the case of condition (16), we have estimated Idv^/dtl with £ 2 ^ |v t I . 

Now, |u 2 | approximates Iv 2 l, and (v 3 | > lu 3 | for most rising or falling 
balloons. Thus, ~ 
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( 27 ) 


n B ^l 1 R B 1V l‘ 

«2 iv 3 i + ifi 3 v 2 i ~ fi 2 iu 3 1 + m 3 u 2 i ’ 

so that condition (26) is satisfied because condition (2i) was found to be valid. 
Thus, it appears that the horizontal Coriolis terms in equation (6) can be 
safely neglected against the horizontal balloon acceleration terms referenced 
to the relative frame of reference and the vertical Coriolis terms can be 
neglected against the Archimedean buoyancy term. 


LINEAR PERTURBATION EQUATIONS 


In the standard linearized balloon response problem, one considers the 
response of the balloon to infinitesimal wind perturbations superimposed on 
an unperturbed wind state which is spatially and temporally invariant. The 
assumed infinitesimal nature of the wind perturbations permits one to assume 
that the balloon responses are infinitesimal and thus allows one to pose the 
balloon response problem in a linear context. The motivation for this pro- 
cedure is the simplicity of the resulting equations which facilitate a relatively 
straightforward calculation of the response properties of the balloon. The 
analysis herein neglects the Coriolis terms in equation (6) . 

We assume that the wind field can be represented in terms of a spatial- 
ly and temporally invariant mean component and a superimposed wind per- 
turbation which can depend on both space and time; 

u. (x^, t) = u\ + u' (x^, t) . (28) 

Henceforth, overbar and prime indicate a mean or unperturbed state and a 
superimposed fluctuation. The assumed mean wind field executes horizontal 
flow, so that 

u 3 = 0 . (29) 

The balloon velocity fluctuations are represented in a similar way; i.e. , 

v. (t) = v. +v|(t) . (30) 
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Unlike the assumed mean wind field, ^ 0. It is assumed that the mean 
horizontal motion of the balloon is in equilibrium with the environmental wind, 
so that 


v. = u., i = i, 2 . (31) 

l r 

The aerodynamic lift coefficient vector is a time -dependent vector. The mean 
aerodynamic lift coefficients are taken to be equal to zero, 



(32) 


Thus, the life coefficient vector is assumed to be a perturbation vector. 

Fichtl, DeMandel, and Krivo [ 7] have shown for a sufficiently fast rising or 
falling balloon that if the horizontal lift coefficients and the wind and balloon 
velocity perturbations are of first-order smallness, then the vertical lift coef- 
ficient is a second-order quantity. In other words, to within first order, the 
lift force acts in the horizontal plane for a sufficiently fast mean rise or fall 
rate. Because the lift coefficients can vary in time, it is reasonable to allow 
for the drag coefficient to depend on time. Thus, we shall assume that 


C 


D 




(33) 


It is assumed that the mean or unperturbed quantities satisfy the equations of 
motion, equation (6) , so that in view of the above assumptions, the horizontal 
equations of motion of the balloon are identically satisfied. The vertical equa- 
tion of motion for the balloon reduces to 



lv 3 I v 3 - g (m - m Q ) = 0 


(34a) 


Substitution of equations (28) , (30) , and (33) into the balloon equations of 
motion [equation (6) ], utilization of conditions ( 29) , (31), (32), and (34a), 
and the neglect of second- and higher-order terms yields 


dv! v! / 3u! 

-« + lir* 

u! 

+ y ( i) TjT + C 2 - y( i)] 



au! > 

i _ *■; 


X 

K > 

) + “ v 3 ax 3 

00 

>_ 

- c* 

.a., A 

TC d L1 

1 TC d 


(34b) 
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where 


y(i) = 



i=i, 2 

i = 3 


2 (m + m ) 
E 

P AC d ^ 3 1 


(m + m ) |v 3 

a 

lm - ml g ’ 


M = 


m+m 
o a 

m + m 

a 


> 


and 

m 

a 

a = 

m + m 

a 


(35) 


(36) 


(37) 


(38) 


The Einstein summation convention does not apply to terms containing y(i) . 

In deriving equation (34b) , we have assumed that all velocity perturbations and 
drag and horizontal lift coefficient fluctuations are first-order quantities, 
while the vertical lift coefficient is a second-order quantity. 


For a sphere, the apparent mass is related to the displaced air mass 
through the relationship 


m 

a 



(39) 


so that the parameters p. and a can be expressed in terms of one parameter, 
e; thus, 


3 

M “ i + 2 e * 


(40) 


a 


i 

i + 2e 


(41) 


where 


e 



(42) 
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Equations (34a) and (34b) are the basic equations of the linear response 
problem. To solve equation (34b) , one must specify the environmental wind 
perturbations. In the special case in which u! is a function of x 3 , only the 

term involving p in equation (34b) vanishes identically. This means that, in 
the context of the linear response problem, the dynamic buoyancy effect 
vanishes when the wind perturbations have only vertical variations. 

FOURIER INTEGRAL REPRESENTATION OF WIND FIELD 
AND AERODYNAMIC COEFFICIENTS 

We hypothesize that the wind perturbations can be represented with a 
Fourier integral as follows: 

u [ ( x j» *) = I f f f \ ( K m > w ) \~ wt ) dKdco , (43) 

J —00 

where ( A ) denotes Fourier amplitude and ck^ckj ck 2 d/c 3 . The quantity i in the 
exponent is the \T^1 and whenever i occurs as a coefficient it is understood 
to be this quantity; otherwise it is an index. The quantity k denotes the k 

component of the radian wave number vector and u> is the radian frequency. 

All quantities which involve u! in equation ( 34b) are interpreted to 
mean that we perform the indicated operation on ul^ and then evaluate at 
x^. = X^(t) . To do this, we must perturb the Lagrangian coordinates of the 
balloon. Integration of equation ( 7) yields 

t 

X^t) = / v k (t) dt , (44) 

where, without loss of generality, we have taken the initial position of the 
balloon to be at the origin of the coordinate system [X^ (0) = 0] . Substitution 
of equation ( 30) into equation ( 44) yields 

x k (t) = \ t+x k(‘> • (45) 

where 

t 

X k (t) = / v i (t) dt ‘ (46) 
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The quantity X^(t) should be a perturbation quantity of first or higher order 

of smallness because the sign of v’ (t)’will fluctuate in time, thus tending to 

cause the positive and negative contributions to the integral to cancel. We 
will assume that X^(t) is a quantity of first-order smallness. 

Substitution of equation (43) into the right side of equation (34b) and 
^valuation at the location of the balloon given by equation ( 45) yield 


dv' v! 

-j-+y<i)Y = E 2 - y( i) ] 



c; . - 6.. 

Li i3 




- //// 

— oo 


i [p(u k K k - co) + a 


V 3 S 1 + 



(47) 


A / \ i(/c, v, - o>)t , , 

x u. (k., w) e k k d« do) , 

1 1 

where we have neglected second- and higher-order terms in fluctuation 
quantities. Let us set 

S2 = K k v k - w , (48) 


so that equation ( 47) can be expressed as 


dv! v| 

~Y +r(i) Y = [2 - y(i)J 






oo 

+ / B. (S 2 ) e lfit di2 

— OO 


(49) 


where 


i (S2) = ^ Iff { 1 [( ^“ a) + 

A ( 

X u. u lf 


^2 — £2 j 


(50) 


K 2, — = , W J d*i cLk 2 dw 

v 3 / 
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and 


The quantities C^. and are functions of time which can also be 
represented as Fourier integrals: 


C 


t 

Li 



dp 


(51) 


and 



C D e iat d C 


(52) 


Substitution of equations (51) and (52) into equation (49) yields 


dv’ v* °° 

k , /M k _ r , _v lfit 

^ y( W ip ~ J ( 12 ) © dp , 

—CO 


(53) 


where 


V V 

<V°> =B k (n) + [2-y(k)J -Jr £ Lk -e k3 -Jr c D 


(54) 


TC 


D 


TC. 


D 


FORCED RESPONSE OF BALLOONS 


To obtain the forced response of a balloon, we assume that 


v L(t)“ / \ (P) e lS2t dp 
-00 


(55) 
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Substitution of equation ( 55) into equation ( 53) yields 


v, (fi) = 


y(k) + ifit G k (fi) 


( 56 ) 


Thus, the forced response of a balloon could be calculated by substituting 
equation ( 56) back into equation ( 55) to yield 


T k« - / 


T 

y(k) + ifit 


G k (fi) e lS2t dfi 


(57) 


The power spectra of the balloon velocities are useful statistics for 
describing the properties of the balloon motions. To obtain the generalized 
power spectrum of v^(t) , we multiply equation (56) by its complex conjugate 

evaluated at fi so that 


<!\< c > % (n '*> 


T 2 <G k (fi) G* (fi')> 

•/(k ) + iy(k) T (fi - fi ') + fi fi f T 2 * 


(58) 


where <( )> denotes the ensemble average operator. Substitution of equation 
(54) into equation (58) yields 


<v fe (fi) v* («')> = 

+ [2 - y( i) ] 



T 2 <B k (fi) B* (fi’J> 

y*(k) + i y(k) T ( fi - fi *) + fi fi ’ T 2 

2 <c Li (n) cj.(fi*)> 

y2(k) + i y(k) T (fi - fi 1 ) + fi fi* T 2 


6 i3 



<^C d (J2) 

y 2 (k) + y(k) T ( fi - fi *) +~S2~TP T 2 


9 


(59) 


where we have assumed that the Fourier amplitudes C T ( fi) and C ( fi) are 

-Lik D 

uncorrelated with B ( fi) . The second and third terms on the right side of 
equation (59) are the contributions to the balloon response spectra from the 
fluctuations in the aerodynamic lift and drag coefficients. These terms have 
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been analyzed in detail for the Jimsphere balloon [5] by Fichtl, DeMandel, 
and Krivo [7] . The first term which we shall denote as S(fi, 12') is the 
contribution to the balloon response sepqtra from the wind perturbations along 
the balloon track. The remainder of this report will be concerned with an 
analysis of this term. 

Substitutuion of equation ( 50) into S ( 12 , 12') yields the balloon response 
spectrum in terms of the spectrum of u| : 

S k (fi, «')=/... f H fe $ ^ k u k 2 , , w, k\, K 2 , ^ ~ 

x d«j dfc 2 dw d#cj dfc 2 do; *, 

where 

-^(k) » [(n - a) toj + an] [(M ~ a) fll + an 'IT* + i y(k) Hu - a) (O t - nj) + a( a - n')]T 2 
k y 2 (k) + iy(k) T(t2 “O'J + fi fi'T 2 



$ 


k 


1 



* 2 » 


12 ~ 12 i 
v 3 



(62) 


The spectrum S k (l2, 12') is the generalized power spectrum as defined 
by Bendat and Piersol [ 8] . Another useful representation of the spectrum is 


¥ (12, t) = f S( 12, 12’) e l( fi " dl2*. 

— OO 


(63) 


This spectrum is the nonstationary spectrum as defined by Kinsman [ 9] in his 
analysis of ocean waves. 

STATIONARY BALLOON RESPONSE 

If the balloon response is statistically stationary, then ^ k must be 
time invariant, which means S( 12 , 12 ') must be of the form 
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S(fi, fi») = ¥ (fl) <5 (a - fi r ). 
k 


( 64 ) 


This implies 


* / n-S2i i r _ / fl-fli \ 

*k \ K i’ *2, “Hgp. "» K i’ *2, “ 8 » w 7 =# kV <1 ' * 2 ’ ~VT^’ U ) 


x 6 (kj - kJ) 6(«2 - Kp 6 ( J2 - R') 6(co) . 

Substitution of equations ( 64) and ( 65) into equation ( 60) yields 


(65) 


err Y (k) + t(M - a) fii + an ] 2 T 2 . / n - Ri \ , , , 

*k (n> -JJJ TlioTTarp M** *»• > “ )<"*< ^ d “- 

' ' ( 66 ) 

This result is the stationary spectrum for the k component of balloon velocity 
resulting from a statistically stationary and homogeneous wind field. 


WIND FIELD WITH ONLY VERTICAL VARIATION 


The result given by equation ( 66) shows that the spectrum of a com- 
ponent of balloon velocity is related to a three -fold integral of the four- 
dimensional spectrum of the corresponding component of the wind. In the 
special case in which the wind only has vertical variation, we must have 

*k (**’ * 2 ’ ~ » ") = , (67) 

so that equation (66) reduces to 


* k (fl) = 


y 2 (k) + (otftT) 2 
y 2 (k) + ( S2T) 2 




( 68 ) 


The quantity 4> ,( ) is the one -dimensional wave -number spectrum of the 

K 

wind along the vertical transformed to a frequency spectrum via the trans- 
formation 


R = v 3 k 3 . (69) 

A detailed account of equation (68) can be found in Reference 6. 
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RESPONSE OF BALLOONS TO ATMOSPHERIC FLOWS IN 
WHICH TAYLOR'S HYPOTHESIS IS VALID 


In many instances, atmospheric turbulence and other types of stochas- 
tic flo.ws seem to obey Taylor* s frozen eddy hypothesis [ 10] , which states that 
local changes are produced by advective effects or 

W = \ K k * (70) 

This means that the spectrum of the wind must be of the form 


$ 


k 



*2 


12 ~ £2 1 
V S 



« 2 » 



(71) 


so that integration of equation ( 66) with this representation yields 


V C, = 


y^(k) + ( <r 12 T) 2 
y z (k) + (12T) Z 


*k (a) 


The quantity <$ ( fi) is a one -dimensional spectrum given by 

OO OO 

*k (S2)= ff *k \ Kl ’ * 2 ’ ^ ) dK i dK2 • 


(72) 


(73) 


The result given by equation ( 72) is identical to the form of the result given 
by equation (66) . The assumptions, however, which lead to the two results 
are different. The result [equation (68)] is a special case of equation (72) 
and was obtained by assuming the wind field can vary only in the vertical. 

This assumption is rather restrictive so that the result given by equation (68) 
is of only a limited interest. On the other hand, the result given by equation 
(72) is intriguing in that, except for the basic assumptions of the perturbation 
scheme, we have only assumed that Taylor’s hypothesis is valid at each level 
(eddies are transported by the mean wind) . Taylor’s hypothesis is not a very 
restrictive assumption and is commonly invoked in studies of atmospheric 
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turbulence. The hypothesis appears to be valid for a wide variety of situations 
[ 10] . The result [equation ( 72) ] says that, if the mean wind is height- 
invariant, then it is possible to compute the one -dimensional spectrum of the 
k component of the turbulence velocity vector directly from the spectrum of 
the corresponding components of balloon velocity. The result offers many 
possibilities with regard to turbulence measurements with balloons. If the 
mean wind is not reasonably constant with height or if the Fourier components 
propagate relative to the mean flow, thus voiding the application of Taylor's 
hypothesis, the simple result given by equation (72) fails. 


Evaluation of equation ( 43) at the location of the balloon and utilization 
of equation (70) yields 


OO 

u! (v 3 t) = f 

— oo 


oo 

// \ U k > d «i d *2 

-00 


e iv 3 *3 t 


d K 3 


(74) 


where we have put 


A 

u. 


(Kj^, 0>) = u j(^) <5 (ctJ - U. K.) 

J J 


(75) 


Now, V 3 t = x 3 , so that equation ( 74) states that the horizontal components of 
turbulence (u[ and u 3 ) along a balloon track in a space -time invariant hori- 
zontal mean flow are lateral velocity fluctuations and the vertical velocity 
perturbations also along the track are longitudinal fluctuations if Taylor's 
hypothesis is valid. This means that 4> j and $ 2 * n die context of equation (72) 
are lateral spectra and <f> 3 is a logitudinal spectrum. 


CONCLUSIONS 


The interpretation of balloon data is not a straightforward procedure. 
This is especially true if the balloon tranverses all or a large portion of a 
flow field Fourier component. If the waves move sufficiently slow relative 
to the balloon rise rate and if the wave lengths of the Fourier components are 
sufficiently long, then the balloon essentially measures a mean flow velocity. 
However, if the waves are sufficiently short and have phase velocities on the 
order of the balloon rise rate, then the interpretation of a balloon wind 
measurement becomes difficult. In this report, we have attempted to define 
a complete set of equations for spherical balloon response studies and have 
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If 


calculated the linear response properties of a spherical balloon to a three- 
dimensional time -varying wind field. We found that, in the case of stationary 
and homogeneous flows, the power spectra of the balloon components of ve- 
locity are related to a three-fold integral of the three-dimensional frequency 
spectrum of the corresponding components of the wind modified by a balloon 
system function [equation (66)] . We found that, in the special case in which 
Taylor’s hypothesis is valid at each level, the above-stated integral can be 
evaluated to yield the result that the power spectra of the balloon components 
of velocity are directly proportional to the one -dimensional spectra of the 
corresponding components of air velocity along the balloon track relative to 
the mean flow. The function of proportionality in this case is the system 
transfer function of the balloon which is a known function. This result is 
identical in form to that obtained by previous investigators (noted in the Intro- 
duction) for the case in which only vertical variations of the wind field are 
permitted. Our analysis showed that the horizontal and vertical balloon ve- 
locity spectra are lateral and longitudinal. These results offer intriguing 
possibilities for turbulence research with balloons because Taylor* s hypo- 
thesis is not an unduly restrictive assumption. This analysis shows that it is 
possible, in principle, to obtain longitudinal and lateral spectra of turbulence, 
provided the underlying assumptions of the analysis are reasonably satisfied. 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama, February 11, 1972 
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